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Abstract. 

Segregation of impurities to grain boundaries plays an important role in both the 
stability and macroscopic behavior of polycrystalline materials. The research objective 
in this work is to better characterize the energetics and length scales involved with the 
process of solute and impurity segregation to grain boundaries. Molecular statics 
simulations are used to calculate the segregation energies for carbon within multiple 
substitutional and interstitial grain boundary sites over a database of 125 symmetric 
tilt grain boundaries in Fe. The simulation results show that there are two energetically 
favorable grain boundary segregation processes: (1) an octahedral C atom in the lattice 
segregating to an interstitial grain boundary site and (2) an octahedral C atom and a 
vacancy in the lattice segregating to a grain boundary substitutional site. In both cases, 
lower segregation energies than appear in the bulk lattice were calculated. Moreover, 
based on segregation energies approaching bulk values, the length scale of interaction 
is larger for interstitial C than for substitutional C in the grain boundary (w 5 A 
compared to w 3 A from center of the grain boundary). A subsequent data reduction 
and statistical representation of this dataset provides critical information such as about 
the mean segregation energy and the associated energy distributions for carbon atoms 
as a function of distance from the grain boundary, which quantitatively informs higher 
scale models with energetics and length scales necessary for capturing the segregation 
behavior of alloying elements and impurities in Fe. The significance of this research 
is the development of a methodology capable of ascertaining segregation energies over 
a wide range of grain boundary character (typical of that observed in polycrystalline 
materials), which herein has been applied to carbon segregation to substitutional and 
interstitial sites in a specific class of grain boundaries in a-Fe. 

Keywords: Impurity segregation, Segregation energy, Grain boundary, Molecular 
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1. Introduction 

The computational design of future alloys will greatly depend on our ability to 
understand and quantify nanoscale phenomena in metallic material systems. For 
instance, impurity segregation to grain boundaries (GBs) in alloys can have a profound 
effect on underlying microstructural processes, which can subsequently be detrimental to 
mechanical properties in polycrystals, e.g., hardness, toughness, and fracture behavior 
[HGJEIIHEJEIITIIEJIHI HSl E]- On the other hand, in some cases, atom segregation 
to GBs can actually be beneficial for macroscale material properties, e.g., by forming 
intermetallics, strengthening GB cohesion, or preventing grain growth [Tij IT2" } [T3"l fH^ [To] . 
Segregation also plays a role in GB decohesion. For instance, Yamaguchi et al. recently 
showed that S segregation to Ni GBs leads to a reduction in GB tensile strength by 
an order of magnitude p2j. Moreover, Solanki et al. found that certain H defects are 
favored at a-Fe GBs and that these species affect the cohesive GB strength jTTj. Since 
the presence of impurities and atoms at GBs can have such an acute impact on many 
material properties, understanding their interaction with and segregation to GBs and 
other lattice defects is crucial to the design of future materials. 

One potential application of work in atomic segregation is nuclear materials. 
Nuclear material design is also dependent upon understanding the segregation of 
impurities and defects within cladding materials. Radiation damage, through cascade 
events, ultimately results in numerous vacancies and interstitial atoms within the lattice. 
Impurities within the material then tend to diffuse with the vacancies or interstitial 
atoms as they attempt to return to equilibrium positions in the lattice [TU EH 120] . 
Such non-equilibrium radiation-induced segregation has a profound effect on material 
properties due to the accelerated segregation kinetics in comparison to the typical 
kinetics in thermal equilibrium [201 |21| 122] . Moreover, since many cladding materials 
are polycrystalline and GBs are significant sinks for defect and impurity segregation, 
understanding impurity segregation to GBs is crucial to nuclear material design. 

A number of studies have experimentally characterized the presence and effect of 
impurities on GBs in various materials pQ|EllI3lE3ll21E5ll23l2ZlEEll22ll3niEIll32ll3^. 
For instance, Lejcek used Auger electron spectroscopy (AES) to show that segregants 
are equally distributed between fracture surfaces in symmetric tilt grain boundaries 
(STGBs) and distributed unevenly for asymmetric boundaries in Fe-Si bicrystals [28J. 
Furthermore, Lejcek et al. comprehensively classified [100] tilt GBs in a-Fe into special, 
vicinal, and general categories using AES measurements of GB segregation [27]. Such 
studies have also proven useful in GB engineering. Recently, Kobayashi et al. used 
electron backscatter diffraction (EBSD) and orientation imaging microscopy (OIM) 
to show that intergranular embrittlement caused by sulfur segregation in nickel can 
be lessened by developing an optimal GB microstructure [29]. Moreover, EBSD 
experiments of Al-Zr alloys have shown that GB sites in immobile twist GBs have 
a much higher degree of segregation than at mobile tilt GBs [23]. Researchers have 
also begun to use high-resolution transmission electron microscopy (TEM) and Local 
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Electron Atom Probes (LEAP) [301 ED E21 [33] to create three-dimensional atom-by- 
atom representations of solute segregation at GBs and characterize their concentrations. 
For example, Taheri et al. utilized a method that combined EBSD and focused ion 
beam milling specimen preparation with LEAP to measure solute segregation at GBs 
in an Al alloy [30J. Furthermore, LEAP has been utilized by Isheim and colleagues 
to illustrate the reduction in impact toughness in low-carbon steels as a result of the 
combined segregation behavior of C, B, S, and P [21] ■ While critical experiments provide 
valuable insight into solute segregation to GBs, techniques that aim to probe how 
atomic structure impacts segregation are often difficult to perform, expensive, and very 
time intensive. Additionally, these sorts of experiments have yet to be used to study 
large numbers of boundaries with varying GB character, typical of real polycrystalline 
materials. 

Modeling and simulation of segregation to GBs at the atomic scale can also 
provide valuable insight into segregation processes in polycrystalline materials JTUJ 

LIHl EHl |36l ETl EHl [391 SQl SH S21 Sal SH SHI SHI SZl SHI SHI IHQl EH [52], [53]. 

Typically, modeling and simulation of GB segregation at the nanoscale uses ab initio 
simulations [121 HH1 IM1 ESI EH1 EHl EH SOI SH S21 S31 or molecular dynamics (MD) 
[SI SSI S6l SH SH S9l ETJI EH [S3]. Ab initio calculations are often used to study 
the electronic effects of solute presence at GBs and their influence on cohesive strength. 
For instance, Liu et al. investigated the preferred site of Mg segregation at Al GBs and 
determined that Mg forms weaker metallic bonds with Al atoms in the GB region and 
decreases the cohesive strength of the GB [37]. Wachowicz and Kiejna |4QJ studied the 
effect of substitutional and interstitial N, B, and O impurities at an Fe GB and found 
that N in both positions and interstitial B are embrittlers while O in both positions and 
substitutional B enhance GB cohesion. The segregation energies and cohesive effects of 
twenty impurities and alloying elements at a Zr twist GB were calculated by Christensen 
et al., who showed that most elements have an adverse effect on GB cohesion, with Cs 
being the most embrittling [38]. These techniques, however, can be computationally 
expensive and have typically been used only for a few GBs. On the other hand, MD 
studies often use empirical or semi-empirical interatomic potentials fit to ab initio and 
experimental properties. These simulations are much less expensive than their ab initio 
counterparts but are limited by the accuracy or availability of interatomic potentials. 
Nonetheless, MD simulations are increasingly being used to study GB segregation in 
both fee and bec materials. Millett et al. investigated the impact of dopants at a 
Cu GB and concluded that, for a particular concentration of each dopant atomic size, 
the thermodynamic driving force for grain growth could be eliminated [45J. Lezzar et 
al. concluded that the driving force for intergranular segregation in Ag(Ni) and Ni(Ag) 
systems can be primarily attributed to the atomic size effect [5D]. While MD has been 
more commonly used for fee materials, such simulations have also provided insight into 
GB segregation in body-centered cubic (bec) Fe as well [SU |52j [53]. For instance, Gao 
et al. used MD simulations to show that, at a-Fe GBs, He binding energy increases 
with excess volume and binds to GBs more strongly in interstitial positions than in 
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substitutional ones [52]. Additionally, Malerba et al. modeled displacement cascades 
in an Fe-Cr system with MD to show that a large percentage of Cr atoms are located 
in interstitial clusters, which may greatly reduce the mobility of interstitial loops when 
compared to pure Fe [53] . 

While MD simulations are much less expensive than ab initio simulations, very few 
simulations consider a large number of GBs in their analysis of GB-related properties. 
GBs have five degrees of freedom associated with them (plus three associated with 
translation at an atomistic level), and many experimental methods have begun to 
measure the GB character in terms of these degrees of freedom [54"] 1331 [56] for GB 
engineering purposes. However, in nanoscale calculations, only a few studies have 
explored fifty or more GBs in their analysis of nanoscale properties. Tschopp and 
McDowell have shown that asymmetric tilt GB systems in Cu and Al facet into the 
structural units of their corresponding symmetric tilt GB counterparts [571 EE1 EH] 
and that the grain boundary structure results in very different dislocation nucleation 
properties and mechanisms [60] EH E2, E3] . Holm et al. calculated energies of 388 GBs 
in Al, Au, Cu, and Ni, and observed that the GB energy scales with the shear modulus 
and that boundaries with significant stacking fault character correlate with the stacking 
fault energy [18]. The classic work of Wolf has shown that, for several Mo and Fe 
GB systems, GB energy correlates nearly linearly with volume expansion per unit area 
[49J. The recent work of Tschopp et al. used >150 Fe STGBs to demonstrate that, 
based on formation energies, self-interstitial atoms display a larger energetic driving 
force for binding to GBs than vacancies do [HJ 163] . Clearly, a similar methodology 
using molecular statics and dynamics simulations that can analyze how segregation in 
a-Fe systems is influenced by GB character would be valuable to understanding GB 
segregation and, perhaps, to engineering materials by increasing beneficial GBs while 
decreasing detrimental GBs. 

In this work, the research objective is to quantify the energetics and length scales 
associated with C segregation to Fe GBs. The Fe-C system is chosen as an ideal 
system because C is known both experimentally E7J EE] and computationally 
[69] 170] 171] 172] 173] to segregate to the boundaries, dislocations, and surfaces; moreover, 
a number of experimental studies have shown that C increases the GB cohesive strength 
in a-Fe [661 EH EH] . The methodology used here provides a means for simulating how 
GB character impacts the segregation of C to a large number of Fe GBs. In this work, we 
utilize an interatomic potential [71] specifically formulated to capture the energetics of C 
interactions with point defects in a-Fe, which is in good agreement with ab initio results 
[73]. This paper is outlined as follows. Section [2] describes the simulation methodology 
used to simulate and calculate segregation data. Section [3] discusses the results of the 
simulations and their significance for modeling GB segregation. Section [4] discusses our 
results, particularly addressing how the present methodology may need to be extended 
to better model the interaction between GBs and interstitial atom species. Section [5j 
summarizes this research and provides conclusions based on our results. 
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2. Simulation Methodology 

In this work, the segregation energy associated with a single C atom was calculated at 
sites within or around a-Fe GBs. While it is well known that C occupies octahedral 
interstitial sites in the perfect single crystal a-Fe lattice, here we examined a few different 
scenarios for C. As a first order approximation, we utilize sites formed on the initial GB 
lattice. This process is meant to mimic the restructuring at the boundary that occurs 
through an interaction with a vacancy and then a subsequent occupation of an available 
site in the restructured boundary (at the exact location of the initial vacancy). The 
hypothesis is that at a GB that has undergone some restructuring due to interactions 
with point defects, such segregation processes and sites may be energetically favorable. 
This may be a reasonable assumption given that DFT calculations have shown that 
another interstitial atom, N, has very similar formation energies in both substitutional 
and interstitial sites within a £5(210) GB [40J. Moreover, in an ultra-low C bake- 
hardening steel sheet, three-dimensional atom probe measurements have found that C 
atom concentrations at the GB can be more than 200 times that in the bulk (68]; hence 
the possibility for C to segregate to both interstitial and substitutional sites is probable. 
In Section |4| this assumption will be further compared with interstitial sites within the 
GB lattice (as opposed to sites directly on top of the GB lattice) for a few GBs, based 
on starting coordinates obtained from a Voronoi tesselation of the simulation cell. The 
first-order process used to calculate the segregation energies of C in a-Fe is as follows: 

(i) A GB is selected from a GB database that contains 125 STGBs (50 (100), 50 (110), 
25 (111)). 

(ii) A GB site (within 15 A) is chosen and a C atom is substituted for the Fe atom at 
this site. 

(iii) A molecular dynamics code (LAMMPS |76j ) is used to minimize the energy of the 
GB with the substitutional C atom. 

(iv) The GB, site position, and calculated segregation energy of the substitutional C 
atom are stored. 

(v) The process is repeated for all sites within 15 A of the GB center and for all GBs 
within the GB database. 

The Hepburn and Ackland Fe-C interatomic potential [74J is used to model the Fe 
GBs and their interaction with the subsitutional C atom. This potential is based on the 
embedded-atom method (EAM) formalism [TTJ [78] and is in agreement with Density 
Functional Theory with respect to the energetics pertaining to interactions between 
C atoms and Fe self-interstitial atoms, vacancies, and other C atoms. Unlike prior 
Fe-C potentials, the Hepburn-Ackland Fe-C potential was the first EAM potential to 
correctly capture covalent bonding of two C atoms within a vacancy. Moreover, previous 
EAM potentials showed strong binding of C to overcoordinated defects, such as self- 
interstitial atoms, whereas the Hepburn-Ackland potential correctly captures the strong 
repulsion between orvercoordinated defects and C, in agreement with ab initio results. 
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This repulsion can be important for interactions between C atoms and GBs. Last, this 
potential has been simulated at temperature yielding dynamics and mechanism^] in 
agreement with ab initio results [75]. For instance, Terentyev et al. [79] recently used 
this potential to investigate the influence of C atoms on the stability and migration of 
small clusters of point defects and found that C atoms have an attractive interaction 
with vacancy clusters containing fewer than four vacancies. This potential provides a 
reasonably accurate representation of the Fe-C system and is deemed appropriate for 
studies of single C atoms within the bcc Fe lattice. 

The segregation energy is calculated for C as a function of position at each site 
within 15 A of the GB. For each GB structure, an Fe atom at a particular site a 
is replaced with a C atom and the simulation cell is relaxed using the Polak-Ribiere 
conjugate gradient energy minimization process. The total energy of the simulation cell 
is calculated and the process is repeated for each atomic site within each GB in the 
database. The segregation energy calculations follow a similar approach to others, e.g., 
Liu et al. [37]. The segregation energy associated with a C atom at site a is calculated 
with 

h seg ~ [^tot ~ &tot J ~ y^bulk ~ h bulk J I 1 ) 

where Ef f' Csub and Ef Q f are the total energies of the GB structure with and without 
the solute substitution. E^ lk Caub and E^ lk are the total energies of a single crystal bulk 
Fe simulation cell with and without the substituted C solute. The bulk energies used a 
10a xl0a xl0a bcc cell with 2000 atoms. Hence, E^ lk is equal to 2000£f e [E F c e = 4.013 
eV, i.e., cohesive energy of Fe) and E^f sub = 2000£f e + 0.391 eV. These bulk energies 
are subtracted in Eq. [l]to remove the effect of substituting the C atom. Hence, Eg" g « 
represents that substituting a C atom into site a in the GB simulation cell results in an 
equivalent energy difference as substituting the C atom into a perfect bcc Fe lattice. As 
with prior work, a negative value of E^" represents that it is energetically favorable for 
C to bind to site a compared to the bulk lattice. It should be noted that the segregation 
energy, E^ g , as defined is equal to (—1) times the binding energy E^ a . Using the same 
terms as in Eq. [TJ the binding energy is typically defined as 

= + Eftf**) - (El?^ + E& h ) = -Eg a g . (2) 

The segregation energy in Eq. [I] will be used for the subsequent analysis of a C atom in 
substitutional GB sites. Later, in Section |4j the segregation energy equation is discussed 
in the context of two different scenarios: (1) a C atom in an octahedral site and a vacancy 
in the bulk lattice combining into C at a GB substitutional site and (2) a C atom at 
an octahedral site occupying GB interstitial sites. That is, the present definition of 
segregation energy in Eq. [I] is modified for the Fe-C system to account for the fact that 

| The migration energy, E m = 0.887 eV, is in good agreement with the ab initio results of Domain 
et al. [75], E m = 0.902 eV. Additional dynamic simulations by Hepburn and Ackland [74] at 1400 K 
with a single C atom in an Fe single crystal lattice show that C exclusively migrates from octahedral 
to octahedral sites through the tetrahedral sites. 
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C does not occupy substitutional sites within the bulk lattice. This energy difference 
is 0.930 eV; hence, E^" g < —0.930 eV is necessary for Scenario 1 (octahedral C and 
vacancy to GB substitutional C) to be energetically favorable. The method outlined in 
this section was used for each site in all 50 (100) STGB, as well as 50 (110) and 25 (111) 
STGBs. For each GB, the segregation energies were calculated as a function of atomic 
location. 

3. Simulation Results 

3.1. Grain Boundary Structure and Energy 

The GB structure database used in the simulations herein contained 50 (100), 50 
(110), and 25 (111) STGBs. Bicrystal simulation cells with three-dimensional periodic 
boundary conditions were used to create the database [EDI ESI EH EH]- To remove any 
possible interaction between the two boundaries, a minimum distance of 12 nm was 
used between them during generation. As with past work [HHl E7J EH] , an atom deletion 
criterion along with multiple initial configurations with various in-plane rigid body 
translations were utilized to accurately obtain optimal minimum energy GB structure 
via the nonlinear conjugate gradient energy minimization process. 

The structures and energies of STGBs may be important to understand the 
interaction between C atoms and the boundary. To examine the range of GB structures 
and energies that might be seen in polycrystalline materials, different grain boundaries 
from several GB tilt systems were used in the present simulations. The database used 
in this work is an expanded version of that first utilized in Tschopp et al. [M]- The 
(100), (110), and (111) STGB systems chosen have several low order coincident site 
lattice (CSL) grain boundaries (e.g., E3, E5, E9, Ell, and E13 boundaries), as well 
as both general high angle boundaries and low angle grain boundaries (6 < 15°). The 
GB energy as a function of misorientation angle for the (100) STGB system is shown in 
Fig. [I] This plot is similar to that found previously in Fe-Cr simulations [81] and similar 
to misorientation-energy relationships found in fee metals [801 [821 [831 EH [85] . The low- 
order CSL grain boundaries for the (100) STGB system (E5 and E13 boundaries) are 
also illustrated in this figure. For the (100) tilt axis, only minor cusps were observed 
in the energy relationship, most noticeably at the E5{310} boundary (990 mJ/m 2 ). In 
addition to many general high angle boundaries, several low angle boundaries (8 < 15°) 
are also plotted. The range of GB energies sampled was 500 mJ/m 2 . 

The GB structure plays an important role on the GB properties [86]. For low 
angle boundaries, the GB is composed of an array of discrete dislocations and the 
corresponding energy can be calculated based on the classic Read-Shockley dislocation 
model. However, at higher misorientation angles, the spacing between dislocations is 
small enough that dislocation cores overlap and dislocations rearrange to minimize 
the energy of the boundary. The resulting GB structures are often characterized by 
structural units [87]. Grain boundaries with certain misorientation angles (and typically 
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Figure 1. (100) symmetric tilt grain boundary energy as a function of misoricntation 
angle [55] • The low-E grain boundaries (E < 13) in each system are identified. 



a low £ value) correspond to "favored" structural units, while all other boundaries are 
characterized by structural units from the two neighboring favored boundaries. An 
example of structural units in the (100) STGB system is shown in Fig. [2j where the 
two £5 boundaries are favored STGBs, and the £29(730) boundary is a combination 
of structural units from the two £5 boundaries. The structural units for the £5(210) 
and £5(310) STGBs are labeled B and C, respectively, in a convention similar to that 
used for face-centered cubic metals [80]. Also, notice that the ratio of structural units in 
the £29 GB can be determined by the crystallographic relationship of the two favored 
boundaries, i.e., £29 (730) = 2 [£5 (210)] + 1 [£5 (310)]. In a similar manner, the two 
£17 boundaries are combinations of the favored B and C structural units and "structural 
units" of the perfect lattice, A and A'. 

3.2. Segregation Energy for (100) Boundaries 

The segregation energies that correspond to the atomic positions in the middle three GB 
structures (Fig. [2]) are shown in Figure [3j AtomEye is used to visualize the simulation 
results [88]. In this graph, the color bar is normalized by subtracting the energy of 
substitutional C in the bulk so that the difference in energy between sites near the GB 
and in the bulk can be easily compared (i.e., atoms colored white have bulk segregation 
energies). For all three GBs, the segregation energy becomes lower as the sites are 
located closer to the GB, meaning that segregation to the GB is favored for substitutional 
C. However, there is not a simple gradient of the segregation energy from the GB center; 
the local structure also plays a pivotal role in the segregation energy. For sites located 
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<100> Symmetric Tilt Grain Boundary Structure 




Figure 2. (100) symmetric tilt grain boundary structures with structural units 
outlined for the £17(410), £5(210), £29(730), £5(310), and £17(530) boundaries [65]. 
Black and white denote atoms on different {100} planes. The different structural units 
are labeled A, B, C, and A'. 

farther from the GB, the segregation energy approaches that of the bulk, as denoted by 
segregation energies close to eV. Interestingly, although the structural units are the 
same between these three grain boundaries, there are some segregation energies in the 
£29(730} that are lower than either of the favored £5(210) and £5(310) STGBs, e.g., 
inside the C structural unit. That is, the elastic interaction between differing structural 
units may produce a different distribution of segregation energies than a boundary 
composed of all the same structural unit. While these trends seem to indicate a driving 
force for the segregation of C atoms from the bulk to the GB, this segregation energy 
needs to be further augmented (E^ g + 0.930 eV) to account for the fact that C lies in 
octahedral interstitial sites in the bulk, as will be further discussed in Section |4| 

Plotting segregation energy against distance from the GB shows information similar 
to that in Figure [3j but provides a convenient method to display the segregation energies 
of the sites in many different GBs at once. The distribution of segregation energies as 
a function of distance for the three GB structures seen in Figure [3] is shown in Figure 
|4| Near the GB, all three GBs show a trend of negative segregation energies at sites 
near the boundary, which is the same behavior reflected in Figure [3} Moreover, notice 
the lack of any segregation energies that are near bulk values within 5 A of the GB 
center for these three boundaries. Figure |4]d is a plot of the same distribution for all 50 
(100) STGBs, which includes both low angle (6 < 15°) and high angle grain boundaries. 
As noted in Figure [4J3, over 10, 000 simulation sites (and atomistic simulations) were 
considered herein. Most of the segregation energies that differ from that of the bulk 
occur between the GB center and about 7 A. While the majority of sites within this 
region have segregation energies less than that of the bulk, there are also a few GB 
sites that have segregation energies that are higher than in the bulk; most of these 
sites tend to be located along the centerline of the boundary. There are a cluster of 
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<100> Symmetric Tilt Grain Boundary Segregation Energies 




£5(310) 229(730) 25(210) 

Figure 3. Segregation energy as a function of site location for substitutional C atom 
in the £5(210), £29(730), and £5(310) boundaries. 

sites around 7-12 A from the GB that have segregation energies lower than the bulk as 
well. There is a subtle difference between low and high angle boundaries. Within 5 A 
of the GB center, low angle grain boundaries tend to have some segregation energies 
that are similar to the bulk values. This is as expected, though. Low angle boundaries 
are composed of dislocations separated by regions of perfect single crystal, which have 
similar segregation energies to bulk energies. 

One way to represent the segregation energies-distance relationship is to bin the 
energies according to their distance from the GB center and to analyze the statistics 
associated with each bin (Figure [5]). Due to the symmetric nature of the GB segregation 
energies as a function of distance (Figure [4]), the absolute value of the distance from the 
GB center was used to provide more data points for the statistical analysis. Furthermore, 
the energies are split into 1 A bins to characterize the distributions and compute 
statistics for sites at a given distance from the GB. An example of the A bin (—0.5 A 
to +0.5 A) is shown in Figure^ along with several statistics: # of boundaries, mean, 
median, standard deviation, and interquartile rangd^j Once the appropriate statistics 
are calculated, a boxplot (Figure^) is used to represent the segregation energy statistics 
in each bin, i.e., the minimum, 25% percentile, median, 75% percentile, and maximum 
segregation energies. In the boxplot, the red line in the box is the median while the top 

§ The interquartile range is denned as the difference between the 25% percentile and 75% percentile. 
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Figure 4. Distribution of segregation energies as a function of distance from the grain 
boundary for (a) the £5(210), £29(730), and £5(310) grain boundaries and (b) all 50 
(100) STGBs. Eg a g < 0.0 eV indicates it is energetically favorable for substitutional 
C to bind to a substitutional site at the grain boundary. However, since C occupies 
octahedral sites in the bulk lattice, E^ g < —0.930 eV is required for an octahedral C 
atom and a vacancy to combine and bind to a GB substitutional site. 



and bottom edges of the blue boxes represent the 25% and 75% quartiles. The whiskers 
extending from the boxes cover the remainder of the range of energies for each bin, and 
the ends of the whiskers denote the maximum and minimum values of the segregation 
energies for each bin. The mean value of the segregation energies in each bin is also 
plotted in green. Boxplots can be very useful for displaying asymmetric distributions. 

The mean segregation energy is lowest with sites close to the GB, as shown in Figure 
[5j and it approaches the normalized bulk value of zero as sites are located farther from 
the boundary. Interestingly, the lowest mean segregation energies actually occur a few 
Angstroms from the center of the boundary. Furthermore, at approximately >8 A, 
the boxes are closely centered about the bulk value, which shows that the overwhelming 
majority of atomic sites display a segregation energy similar to the bulk value. However, 
it is noticed that there are a number of sites with segregation energies significantly below 
the bulk value that still persist up to approximately 11 A. This trend indicates that it 
may be energetically favorable for substitutional C to segregate to sites within 11 A 
of the GB, albeit there is a much larger driving force with decreasing distance from 
the boundary. Additionally, the majority of bins display energy distributions that are 
skewed, usually in the direction of negative energy, i.e., the median is closer to the 
lower edge of the box (mainly for distances less than 7 A). While the median fluctuates 
somewhat, the mean segregation energies - which track with the median - follow a much 
smoother relationship with distance. 
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(a) (b) 

Figure 5. (a) The distribution of segregation energies within 0.5 A of the grain 
boundary center and the associated statistics, (b) Boxplots of segregation energy as 
a function of distance from the grain boundary for all 50 (100) STGBs. The data is 
divided into 1 A bins, and a boxplot is made for each bin. The red lines are medians, 
the blue box ends are the first and third quartiles, and the black whisker ends are 
minimum and maximum values. The mean segregation energy is also plotted in green. 



3.3. Segregation Energy for (110) and (111) Boundaries 

The same process used for the analysis of (100) data in Figures [3]{5] has been repeated for 
the data of (110) and (111) STGB simulations. The distribution of segregation energies 
as a function of distance from the GB for all 50 (110) and 25 (111) STGBs is shown in 
Figure [6] This distribution is similar to that of the (100) STGBs shown in Figure [3Jd. 
However, the minimum segregation energies are much lower than that of (100) STGBs 
and there are fewer sites with segregation energy higher than that of the bulk. 

A statistical representation of the data in Figure [6] is shown in Figure [7j Similar to 
Fig. [HJd, the data has been binned into 1 A bins and the median, quartiles, minimum, and 
maximum values of the segregation energies contained within each bin are shown within 
the boxplots. The mean segregation energy plots trend similarly to that in Figure [5j 
though they display initially lower values close to the GB. The minimum energy whiskers 
again show favorable C segregation sites in most bins: up to 9 A for (110) STGBs and 
up to 11 A for (111) STGBs. 

3.4- Statistical Characterization of Segregation Energies 

Ideally, it would be advantageous to be able to analytically describe the evolution of the 
segregation energies as a function of distance from the GB. Figure [8] provides further 
statistical data for the binned distributions of segregation energies at given distances 
from the GB. The mean and standard deviation of the segregation energy distributions 
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Figure 6. Distribution of segregation energies as a function of distance from the 
grain boundary for 50 (110) and 25 (111) STGBs. Most grain boundary sites within 
8 A have negative segregation energies that decrease with decreasing distance to the 
grain boundary center. 



are plotted in Fig. [8^i. For each of the GB systems, the mean values trend similarly 
between 5 and 15 A, with segregation energy decreasing with increasing distance from 
the GB. However, within 5 A, the mean segregation energy for the (100) STGB system 
(E seg = —0.33 eV) is significantly higher in magnitude than the (110) and (111) STGB 
systems (E seg = —0.63 eV and E seg = —0.74 eV, respectively). Also plotted in Figure 
is the standard deviation of the distributions, which steadily decreases toward zero 
as distance from the GB increases. The decrease is primarily due to the increasing 
number of sites with bulk energy values at distances far from the boundary. For 
normal distributions, the mean and standard deviation would be appropriate statistical 
descriptors to capture the segregation energies. However, the boxplots in Figs. [5] and [7] 
clearly show that the distributions are asymmetric to some degree and may have some 
extreme values or outliers. 

To quantify the asymmetric distributions, the kurtosis and skewness of the 
distributions are plotted in Figure [8Jd. The kurtosis is a measure of how heavily the 
variance of the distribution is affected by extreme deviations, or outliers. Skewness is a 
measure of the asymmetry of a distribution and denotes in what direction a distribution 
possesses a longer tail of values. Figure [8]d shows that the kurtosis is relatively low for 
most bins but becomes very large for some GB systems at approximately 8 to 11 A. This 
is due to most of the segregation energies approaching bulk values (as viewed by the 
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Figure 7. Boxplots of segregation energy as a function of distance from the grain 
boundary for (a) (110) and (b) (111) STGBs. As in Fig. [5} the data is divided into 1 
A bins and a boxplot is made for each bin. The red lines are medians, the blue box 
ends arc the first and third quartiles, and the black whisker ends are minimum and 
maximum values. The mean segregation energy is plotted in green. 



small box, or interquartile range) except for a few negative extreme values. Interestingly, 
in this range, increasing kurtosis correlates with decreasing skewness, which is negative 
for all but a few bins over all three GB systems. The skewness indicates that the 
majority of segregation energy distributions possess longer tails of negative energies; 
the kurtosis indicates when these tails are typically the result of extreme deviations. 
Since the magnitudes of these measures become very large at distances far from the 
GB, a great majority of these sites have the bulk value of segregation energy. These 
four statistical parameters (mean, standard deviation, skewness, and kurtosis) can be 
used to better approximate asymmetric segregation energy distributions, e.g., using the 
Pearson system of distributions. 

4. Discussion 

Carbon was inserted into prior Fe sites at the GB to examine the influence of a wide 
range of GB structures on the distribution and magnitude of the segregation energies 
within the GB region. However, in the perfect single crystal a-Fe lattice, it is known that 
C occupies octahedral sites. Hence, to accommodate the difference in energy between 
an octahedral C atom (and vacancy) in the bulk lattice and this same C atom occupying 
a substitutional site at the GB, the segregation energy associated for substitutional site 
a can be modified, i.e., 

J?C a _ / j? GB > C sub _ TpGB\ _ ( rpFe,C ct _ rpFe \ , ( rpFe,vac _ rpFe \ (o\ 
^seg* — y^tot ^tot J y^bulk ^bulk J "r [^bulk ^bulk J \°) 
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Figure 8. Statistical data for binned segregation energies from Figures 5 and 7: (a) 
mean and standard deviation and (b) kurtosis and skewness. 



where the Fe cohesive energy E^ e is added to the lefthand term in parenthesis to account 
for using the bulk energy from C in an octahedral site E^' k oct in the righthand term. 
Since the difference between Eq. [2]and Eq. [3]is constant, this can be calculated. In Eq. [2j 
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(E^"-E^ k ) = 0.391 eV. In Eq.§ (E^-E^ lk ) = -6.273 eV (solvation energy) 
and (E^ ac - E[ u e lk ) = 5.734 eV. The difference between Eq.gand Eq.gis 0.930 eV. In 
other words, this energy (0.930 eV) must be added to the preceding analysis to compare 
the energetic favorability of C in a substitutional site with that of an octahedral site in 
the bulk lattice. Interestingly, there are a significant number of GB substitutional sites 
that represent a lower energy configuration for a vacancy and interstitial C in the bulk 
lattice (E seg < 0.930 eV). Specifically, the interaction length scale where this process is 
energetically favorable is approximately 3 A from the GB center (total width of ~ 6 A) 
with the largest percentage of favorable sites occurring at the GB center: 7.2% (within 
0.5 A), 5.1% (0.5-1.5 A), 1.1% (1.5-2.5 A), 0.0% (2.5-3.5 A). Additionally, the statistical 
representation of the segregation energy distributions can be used to rapidly quantify 
the probability of lower energy sites within the GB region as a function of distance from 
the GB plane. 

4-1- Substitutional versus Interstitial Carbon in Grain Boundary Region 

Interstitial sites at the GB can be assessed using a similar methodology. There are an 
infinite number of potential interstitial starting positions that could be chosen since there 
is not a set lattice for GB interstitial sites as with substitutional sites. The methodology 
chosen for selecting interstitial sites was based on the Voronoi tesselation method. The 
atom positions within each GB simulation cell was used in tandem with a Voronoi 
tesselation of the cell to generate a list of potential starting positions for interstitial 
sites. In a three-dimensional space, a Voronoi tesselation divides the space into a set of 
space-filling polyhedra that have the following properties: (1) any point located within 
a polyhedron is closest to only one atom, (2) any point on a polyhedron face is equal 
distance to two atoms, (3) any point on a line connecting two polyhedron faces is equal 
distance to three atoms, and (4) any point located on a polyhedron vertice is equal 
distance to four atoms. We have chosen to use the vertices of the Voronoi tesselation to 
populate the set of potential interstitial sites. In the perfect bcc lattice, the polyhedron 
is a truncated octahedron with 14 faces (8 regular hexagonal and 6 square), 36 edges, and 
24 vertices; the Voronoi vertices are located at interstitial tetrahedral sites (i.e., equal 
distance to four Fe atoms). At the GB, however, the polyhedron takes on different 
shapes and the Voronoi vertices constitute sites that are equal distance to four atoms, 
which could potentially be located at the center of GB free volume regions. While a 
perturbation of this technique could be used to identify interstitial octahedral sites in 
the bulk latticdjj] the present technique is deemed sufficient for identifying potential 
interstitial sites in the GB. 

The following Voronoi-based methodology was applied to 50 (100) STGBs (> 

|| For instance, some perturbations of the present technique might be to use the midpoints of the 
polyhedra lines connecting the faces or use the centers of polyhedra faces. In particular, the center 
of the {100} faces formed by the four tetrahedral sites would be in exactly the minimum energy 
octahedral site (0,0.5ao,0.5ao). However, this technique would also include the center of the {111} 
faces (0.25ao,0. 25ao, 0.25ao), which turns out to be a high energy interstitial position for this potential. 
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60,000 sites). First, the distance that each C was displaced during minimization was 
calculated to examine how far each interstitial C moved from its initial site placement. 
This analysis showed that most C atoms (94.5%) were displaced < 0.2 A during the 
energy minimization technique, indicating that C initially placed in the tetrahedral 
sites tends to find a local minimum in energy and does not move to a neighboring 
octahedral site. Even within the GB region, where most (> 99.7%) C atoms displaced 
greater than 0.2 A lie, the maximum distance that the C atom was displaced from the 
initial site was only 0.59 A (< 0.25ao ~~ the minimum distance from a tetrahedral to 
an octahedral site). Hence, the initial positions for the C interstitials identifies local 
minimum energy configurations centered around the tetrahedral interstitial sites in the 
bulk lattice, and the greatest displacements occur within the GB region, as would be 
expected. Again, other perturbations of locations based on a Voronoi tesselation of 
the simulation cell may result in finding interstitial sites with even higher segregation 
energies at the boundaries, but it is anticipated that the present analysis will capture 
the relative influence of interstitial sites segregating to the boundary. 

The results of inserting C atoms at interstitial sites was then analyzed in a similar 
manner to C placed at the substitutional sites. In this analysis, Equation [T] is modified 
such that the segregation energy associated with a C atom at site a, Eseg* is calculated 
by 

E°sk = (^f ^ - Eg) - - E% k ) (4) 

GB C a 

where E tot ' int is the total energy of the GB structure with an interstitial C atom. 
Figure [9|a-d) corresponds to Figs. H and |5j except that the calculated segregation 
energies are for interstitial C using Eq. |4} There are several minimum energy states 
for interstitial C far away from the GB (Fig. [9^a)). The energy of ~ 0.887 eV 
corresponds to the tetrahedral site. However, the present Fe-C potential also has 
minimum energy interstitial sites at other locations, and the Voronoi vertex technique 
did not locate the octahedral site for the C atom upon energy minimization (i.e., 
Eseg 1 = is noticeably absent at large distances). Additional simulations varying the 
location of interstitial positions in a lOaoxlOaoxlOao bcc cell with 2000 atoms show 
that one of the high energy sites is directly between two Fe atoms along the (111) 
direction (0.25a ,0.25a ,0.25a ). In Fig. |9](a), the distribution of segregation energies as 
a function of distance is shown for the same three GB structures as in Fig. ga). The 
GB region again shows both energetically favorable and unfavorable sites, with several 
energetically-favorable interstitial sites having E S eg l of up to —0.5 eV (i.e., a binding 
energv of approximately 0.5 eV). This same trend is also evident for all 50 (100) STGBs 
b)). In fact, based on this plot, the interaction length scale [Eseg 1 < 0) of 
this GB system is on the order of 10 A or less. Binning the data from this plot, the 
distribution of segregation energies within 0.5 A from the GB center (Fig.J^c)) show an 
asymmetric multimodal character with a few peaks centered about the several minimum 
energy interstitial states observed at large distances from the GBs. Approximately 6.9% 
of interstitial sites sampled showed a lower energy than C at an octahedral site in the 
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bulk lattice. The boxplots of interstitial C shows a decrease in the mean segregation 
energy, the interquartile range, and the minimum segregation energy starting at w 5 
A from the GB center (Fig. [9|d)). Interestingly, there is a larger percentage (11.8%) 
of energetically-favorable sites in the second bin (0.5-1.5 A from GB center) than in 
the first bin (within 0.5 A) and this percentage decreases with increasing distance from 
the boundary: 6.9% (0.5 A), 11.8% (0.5-1.5 A), 6.2% (1.5-2.5 A), 3.0% (2.5-3.5 A), 
0.3% (3.5-4.5 A), and 0.2% (4.5-5.5 A). There is a high degree of anisotropy in the 
segregation energies in each bin due to the GB character. Moreover, in contrast to the 
energetic length scales for point defects in a-Fe using the same interatomic potential, 
the calculated length scales of interaction between C and the GB are much lower. 

4-2. Methodology application and extensions 

There are a number of studies that examine the influence of various alloying elements and 
impurities with grain boundaries. For the case of steels, while the interaction between 
GB structures and C represents one example of segregation that has been experimentally 
observed [SHEZIES], many studies also focus on which elements segregate to the GB and 
how these elements may interact with other alloying elements. For instance, elements 
that can impact the properties of steels include, for example, phosphorus, sulfur, 
hydrogen, nitrogen, manganese, silicon, molybdenum, nickel, chromium, antimony, 
and tin [e.g., |89j [90j |9TJ I92J [93] . Also, experiments have shown that there may be 
competitive processes between elements (e.g., C and phosphorus [66]) at the GB that 
will also depend on processing conditions, such as equilibration temperature and quench 
rate. Of course, this is just one example of the complexity of segregation in steels 
due to the numerous interactions of different elements with grain boundaries and the 
kinetics of those processes. The present work may be modified to examine how binding 
behavior between two different atom species (e.g., C-P, C- vacancy) or different atom 
configurations (e.g., H vs. H 2 [17]) are affected by proximity to and within grain 
boundaries. Understanding these interactions may aid in understanding the complexity 
of segregation and its subsequent impact on properties in a-Fe and alloy steels. 

This methodology makes possible the statistical representation of impurity 
segregation to grain boundaries while accounting for differences between GB structures. 
Thus, the results from this method could be used as inputs in other simulations, such 
as the kinetic Monte Carlo technique, mesoscale models, or analytical models. In this 
sense, the information being passed is not just scalar values but is distributions of 
values, which can be used to analyze sensitivity and incorporate variability due to GB 
structure within multiscale models. Clearly, the present application of calculating these 
distributions for segregation of C to grain boundaries represents one such example of 
this concept. Moreover, such studies could be expanded upon to include the the effects 
of temperature with different solutes and solute concentrations. For example, Rittner 
and Seidman [HI] conducted such a study for 21 (110) GBs to calculate segregation 
free energies, entropies, and internal energies for a Ni-Pd system. Their work found 
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Figure 9. Distribution of segregation energies for interstitial sites as a function of 
distance from the grain boundary for (a) the £5(210), £29(730), and £5(310) grain 
boundaries and (b) all 50 (100) STGBs, as in Fig. [4] (c) The distribution of segregation 
energies within 0.5 A of the grain boundary center and the associated statistics, (d) 
Boxplots of segregation energy as a function of distance from the grain boundary for 
all 50 (100) STGBs, divided into 1 A bins and colored as in Fig. [HJ 



a linear relation of segregation internal energies and entropies, which suggests the 
possibility for estimating segregation free energies from internal energies, an easier 
quantity to calculate. Simulations at temperature may lead to an even better prediction 
of segregation behavior to grain boundaries in polycrystalline materials. 
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5. Conclusion 

In this work, we have used molecular statics simulations to investigate the segregation 
energy of a single C atom to thousands of substitutional and interstitial atomic sites 
in 50 (100), 50 (110), and 25 (111) STGBs. A large number of boundaries, including 
general low and high angle GBs, were used in order to account for the variability in GB 
degrees of freedom observed in experimental polycrystalline materials. We can draw the 
following conclusions based upon our results: 

(i) A methodology for calculating and analyzing the segregation energies of thousands 
of sites within a large number of grain boundaries with molecular statics simulations 
has been developed. This method samples different boundaries from a grain 
boundary database and calculates the segregation energy for every grain boundary 
site to acquire segregation statistics. As a first example, C segregation to a- 
Fe boundaries was examined. Both substitutional sites were sampled as well as 
potential interstitial sites, using the polyhedra vertices calculated using a Voronoi 
tesselation of the three-dimensional a-Fe coordinates. Such a methodology is 
warranted given that we found a large degree of anisotropy in sites and segregation 
energies due to varying grain boundary structure. 

(ii) The local structure within the grain boundary affects the segregation energy. As 
an example, the (100) symmetric tilt system is shown where the two £5 grain 
boundaries are both cusps in the energy relationship (Fig. [T]) and contain the 
favored structural units of this system (Fig. [2]) . However, boundaries of intermediate 
misorientations (e.g., the £29 boundary) - which contain combinations of the same 
structural units - do not necessarily have the same segregation energy distributions 
as the £5 boundaries (much lower, see Fig. |3|. 

(iii) For the substitutional C atom case, we found that it is energetically favorable 
for interstitial octahedral C and a vacancy in the lattice to combine within the 
grain boundary at a substitutional site. While this process is highly unfavorable in 
the lattice, there is a region that extends ~ 3 A from the grain boundary center 
where there are favorable substitutional sites for the grain boundaries sampled. The 
largest percentage of favorable sites are directly at the grain boundary center. 

(iv) For the interstitial C case, we found that it is energetically favorable for a C atom at 
an octahedral site in the lattice to segregate to the grain boundary with a maximum 
binding energy of ~ 0.5 eV. The interaction length scale of the grain boundary with 
octahedral C is ~ 5 A from the grain boundary center with the largest percentage 
of favorable sites located within the bin just outside of the grain boundary center 
(0.5 A to 1.5 A from the grain boundary center). 

(v) To quantify the segregation energy distributions as a function of distance from 
the grain boundary, the energies were separated into 1 A bins and characterized 
using several statistical descriptors: quartile values, median, mean, and extreme 
values (see Figs. |5] and [7]). The grain boundary atomic sites have asymmetric 
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distributions of segregation energy with some extreme values that extend over 
10 A from the grain boundary. Furthermore, close to the grain boundary, the 
majority of these distributions are negatively skewed, indicating longer tails of 
negative segregation energies. An analytical model informed by these calculations 
whereby the segregation energy distribution as a function of distance is captured 
using four statistical parameters (mean, standard deviation, kurtosis, skewness - 
see Fig. [8]) is hypothesized for upscaling to higher scale models, i.e., parameters 
necessary for a Pearson system of distributions. 

The significance of this research is not just the calculations of the energetics of C 
segregation in a specific class of grain boundaries in a-Fe, but also the development of 
a methodology capable of ascertaining segregation energies over a wide range of grain 
boundary character typical of that observed in poly crystalline materials. 
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